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Abstract 

We present a bouncing cosmology which evolves from the contracting to the expanding phase in a 
smooth way, without developing instabilities or pathologies and remaining in the regime of validity of 
4d effective field theory. A nearly scale invariant spectrum of perturbations is generated during the 
contracting phase by an isocurvature scalar with a negative exponential potential and then converted 
to adiabatic. The model predicts a slightly blue spectrum, ns > 1, no observable gravitational waves 
and a high (but model dependent) level of non-Gaussianities with local shape. The model represents 
an explicit and predictive alternative to inflation, although, at present, it is clearly less compelling. 



1 Introduction 



Bouncing cosmologies, in which the present era of expansion is preceeded by a contracting phase, 
have been studied as potential alternatives to inflation in solving the problems of standard FRW 
cosmology. The cosmological history is continued to the infinite past, so that all the problems related 
to imposing extremely finely tuned initial conditions in the high curvature regime are avoided. The 
pre Big-Bang scenario pQ and the ekpyrotic/cyclic models are different incarnations of these ideas 

So far, two main problems have prevented these models from becoming serious contenders to 
inflation. 

1. No completely explicit model of a bouncing phase has been found. Both in the pre Big-Bang 
scenario and in the ekpyrotic/cyclic one the bouncing phase lies outside the regime of validity 
of the effective field theory, so that the issue requires some input from a more fundamental 
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theory as string theory. The study of time-dependent solutions in string theory and of the way 
singularities are resolved is not sufficiently developed to describe a bouncing phase or even to 
assess if a bounce is possible altogether [H IH El El l9l [TO]. 

2. No compelling mechanism for producing an approximately scale-invariant spectrum of adia- 
batic perturbations has been found. The pre Big-Bang model predicted a dominant isocur- 
vature mode, which is nowadays vastly incompatible with data. Although some alternatives 
have been studied, scale invariance does not come out naturally [lj. In the ekpyrotic/cyclic 
case it was argued that the same field leading to the fast contraction towards the bounce gives 
rise, with its quantum fluctuations, to an adiabatic, approximately scale invariant spectrum. 
The issue is not completely straightforward because, as discussed, the bounce is not under 
full control so that the way perturbations evolve through it cannot be made explicit. The 
presence of a dynamical attractor simplifies the problem and predictions can be made under 
very general assumptions about the bounce [11]. Unfortunately the result is not scale-invariant 
[12], [131 CD] • O ne cannot rule out the possibility that the bouncing phase has very peculiar 
features which change the result [15], but in this case all predictions will depend on the details 
of the unknown phase. 

Recently there has been progress on both these issues and the purpose of this paper is to put 
these new ingredients together to obtain an explicit and controllable bouncing model with an ap- 
proximately scale invariant spectrum of density perturbations. Let us discuss how to address the 
two problems discussed above. 

1. In a recent paper [16] we studied explicit 4d models of a bouncing phase of the Universe. 
To induce the bounce, the stress-energy tensor must violate the null energy condition. In a 
spatially flat FRW metric this condition corresponds to the inequality H < 0, which is clearly 
violated at the point where we reverse from contraction to expansion: H = and H > 0. 
The violation of the null energy condition is usually associated with catastrophic instabilities 
which make the theory pathological. However this is not true in complete generality, as it can 
be shown with the quite general low energy effective field theory approach developed in [16] . 
An explicit realization can be obtained starting from the ghost condensate model [17J that we 
are now going to briefly describe. In an expanding Universe one expects that a scalar field 
4>{t) evolving in time is progressively driven to rest by the Hubble friction. However, with a 
generic Lagrangian invariant under the shift symmetry (j) — » (j) + const. 

£ = P{{d<j>f) , (1) 

it is easy to realize that another possibility is that (j) goes to a constant, <fi — ► Mg C , in corre- 
spondence of a minimum of the function P. The peculiarity of this solution is that, although 
the scalar keeps on evolving at constant speed, the stress energy tensor is the one of a cosmo- 
logical constant, so that the background metric can be de Sitter or Minkowski. Fluctuations 
around this background solution are healthy and this theory has been studied as a consistent 
modification of gravity in the IR. To understand the relevance of this model in our context 
one can study the stress-energy tensor of the perturbations around the background solution: 
4>(x, t) = Mg C (t + 7r(x, t)). T^v starts linearly in tt with a term [T7] 

T^v = Mg C 7r S^SuO • (2) 
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This shows that a fluctuation with tt < has negative energy! This property opens up the 
possibility of violating the null energy condition. In particular an explicit model of a bounce 
can be obtained adding to the Lagrangain above a suitable potential V ((/)). We will see for 
example that a bouncing phase with a constant positive H can be obtained simply with a 
parabolic potential. In |16j we checked that, with a suitable choice of the model parameters, 
no dangerous instabilities are present. 

Besides giving an existence proof of a smooth bounce, our results allow to study the evolution of 
perturbations across this phase. Nothing exotic happens in these models, so that the general 
conclusions of [TT] apply. Therefore the simple mechanism for producing a scale invariant 
spectrum of perturbations proposed in the ekpyrotic/cyclic scenarios does not work. The 
fluctuations in the field leading the contracting phase have a very blue spectrum, which is 
completely negligible on the scales of cosmological interest. This motivates us to look for an 
additional source of perturbations. 

2. Recently Lehners, McFadden, Steinhardt and Turok [18] pointed out that a second scalar field 
whose energy density is negligible and that evolves along a negative exponential potential 
could be the source of a scale-invariant spectrum of perturbations. These fluctuations can be 
easily converted to adiabatic and therefore match observations. Early studies of isocurvature 
perturbations in bouncing scenarios can be found in [191 120j. 

A possible emerging picture is therefore as follows. The Universe is contracting and its energy 
density is dominated by a field <j>, with a Lagrangian similar to the ghost condensate theory, with 
the addition of a suitable potential. As in the ekpyrotic/cyclic scenario the contraction satisfies 
H -C —H 2 . At a certain point H flips sign until H gets to zero and the Universe starts expanding 
and smoothly connects, after reheating, to a standard FRW cosmology. During the contracting phase 
a second field tp evolves along a negative exponential potential. Perturbations in this "isocurvature" 
field are then converted into adiabatic ones: as eft perturbations are very suppressed, this second 
source gives the leading contribution to cosmological inhomogeneities. 

We are going to study the model starting in Section [2] from the field ifj, which is the source of 
perturbations and therefore the sector of the model giving all the predictions. We will see that, 
neglecting gravity, an exactly scale-invariant spectrum of perturbations is obtained for any negative 
exponential potential. Gravity slightly modifies the spectrum and the resulting tilt is blue. The 
deviation from scale invariance may be very small, but it is not possible to get a red spectrum. This 
prediction is disfavored by recent data [21] and the model will therefore be ruled out in the very 
near future if the detection of a red tilt is confirmed. Additional predictions are the presence of a 
quite high (but model dependent) level of non-Gaussianity (with a local shape dependence) and an 
unobservably small level of primordial gravitational waves (as in the ekpyrotic/cyclic scenarios). In 
Section [3] we study the unperturbed evolution led by the <j> field. It is important to stress that all 
the details are very model dependent and that even the use of a ghost-condensate like theory should 
be considered simply as a working example of an healthy bouncing cosmology. As all predictions do 
not depend on how the contracting and bouncing phases are realized, one can envisage other very 
different possibilities. Some remarks about the exponential form of the potential, which is the root 
of scale invariance in our model, is made in Section H] before the conclusions in Section [5j In the 
Appendix [A] we study the full Lagrangian of perturbations of </> and tp, including mixing terms. 
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2 Exponential is good 



Let us for the moment completely neglect gravity. Consider a scalar ip with standard kinetic term 
and potential V(ip). The unperturbed space independent solution ip(t) will satisfy 

$(t) + V(V(t)) = 0. (3) 

Perturbations around it will be denoted by 5ip(x,t). We want to answer the following question: 
which form of the function V(ip) gives rise to a scale invariant spectrum of perturbations Sip! 
The equation of motion for Sip is 

54>+ [k 2 + V"(^j{t))]5ip = . (4) 

If V" < 0, when the gradient term becomes negligible with respect to the mass term the mode 
"freezes", i.e. the solutions of the equation above are no longer oscillatory. To have a power-law 
solution for Sip(t) in this regime one requires V"(ip(t)) oc t~ 2 . For t < the mass term increases (in 
modulus) and more and more modes freeze. 

As we are assuming V"(ip(t)) oc t~ 2 , equation @ has solutions of the form 

ifj k ~ -J= • F(kt) (5) 

in the conventional normalization. 

Therefore to get a scale invariant spectrum the evolution after freezing must be of the form 
Sip oc t . This time evolution is solution only if 

v"m)) = -f ■ (6) 

With this choice, besides the solution Sip ~ t , we have a decaying solution Sip ~ t 2 . The explicit 
solution with the standard normalization is given by 

<% = -^e~ ikt (l-f). (7) 

Clearly the result of eq. ([6]) could have been anticipated, as eq. (JH) is now the same as the one 
describing a massless scalar living in de-Sitter space (replacing t with the conformal time rj). 

Now we should look at eq. Q as an equation for the function V{ip). Taking the time derivative 
of equation ([3]) we get 

'i> - ^ip = . (8) 

Neglecting a decaying term this gives ip(t) = —2M/t, with M an integration constant with dimension 
of a mass. Thus we have 

ip(t) = -2M\og(-t) + c. (9) 

We can now get rid of time in the second derivative of the potential V"(ip(t)) = —j? and then 
integrate to finally get 

V{ip) = -V e^ /M . (10) 

We thus reach the conclusion that, in the absence of gravity, a field moving along any negative 
exponential potential generates a scale invariant spectrum of perturbations [22]. Notice that while 
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generating say 60 e- folds of scale invariant spectrum the potential varies by e 60 . This implies that 
we really need an exponential potential, and that our result is not an artifact of requiring exact scale 



mvanance 1 ' 



We now want to see whether this way of generating a flat spectrum survives when we put gravity 
back into the game. If the field ip is cosmologically relevant, i.e. if its contribution to the total 
stress-energy tensor is significant, then the picture drastically changes due to its mixing with gravity 
Therefore we are going to assume on the contrary that the energy of ip is negligible so that we 
can forget about its mixing with gravity 0. 

Let us assume that ip lives in a contracting Universe with a oc \t\ p , H = j, with t < and 
p < 1/3. When this last inequality is satisfied, the background solution blue-shifts so fast that gets 
rid of initial curvature, inhomogeneities and anisotropies. In this background the equation of motion 
for ip(t) is given by 

t M K ' 

It is straightforward to check that we still have a power-law solution with 

2M 2M 2 (l-3p) 
T i 2 

The total energy of ip is obviously no longer conserved and it is given by P4> = QpM 2 t~ 2 . This 
contribution will be negligible with respect to the total energy density of the Universe p = 3MpH 2 = 
3p 2 Mpt~ 2 , for sufficiently small M. The solution described above is tuned, as the divergence of ip 
and H both happen at t = 0. This case is clearly not generic, but if we shift the background solution 
H = p(t + to) -1 ) w e see that eq. (fl~2l) is still an approximate solution of (TlTI) going sufficiently back 
into the past: |t| 3> \to\. Let us study the behaviour of perturbations, that are now going to deviate 
from an exactly scale invariant spectrum. 

Tilt of the spectrum. In the presence of gravity, the equation describing the perturbations is 



iP = -— V= 72-^- ( 12 ) 



Sip + y<5^ + 



k 2 2(1 - 3p) 



Sip = . (13) 



t 2 P t 2 

The equation can be simplified using the conformal time r/ and the variable u = aSip. At first order 
in p we obtain 

d 2 u . o 2 — 3p. , 

The solution of this equation with the correct limit at short distance is given by 

u k (v) = ^e^ 1+ ^^^H kv ^- p (-kn) . (15) 

Taking the limit rj — > we get that the function goes, at first order in p, as k~^l 2+p . The tilt is thus 
blue and given by 

n s - 1 = 2p . (16) 



1 Actually we do not know if the perturbations are scale invariant over 60 e-folds, but only in a much 
narrower range. For sure we need an exponential potential in the observable window of ~ 10 e-folds. 

2 The precise regime in which this approximation is justified depends not only on the relative contribution to 
the energy density, but also on the nature of the field that drives the expansion of the universe. In Appendix 
El we shall determine precisely the conditions for our approximation to be valid; for the moment we just 
assume that we are in such a regime. 
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As there is no natural, smooth transition between the rapidly contracting phase when perturbations 
are produced and the bounce, we cannot establish a "natural" value for n s — 1; the deviation from 
scale invariance can be arbitrarily tiny. On the contrary it is quite natural in inflationary models to 
have deviations from scale invariance (of both red and blue kind) of order n s — 1 ~ ± -/V" 1 ~ ± few %, 
where N e is the number of e-folds to the end of inflation. Moreover here we do not have any analogue 
of the slow-roll expansion, so that there is no reason why the tilt should be approximately a constant 
or equivalently why the parameter p should not change in tim^§: in particular the tight constraints 
that we have on p on large scales cannot be extrapolated to short scales. The only sharp prediction 
is that the tilt cannot be red. 

Fluctuations in the ip field will be converted into adiabatic perturbations, which can be de- 
scribed by the usual variable (. We postpone to the next Section the discussion about possible 
conversion mechanisms. The conversion will happen when modes of cosmological interest are frozen 
and gradients are irrelevant, so that the tilt of the spctrum cannot change. At linear order 

where M c is a parameter with dimension of a mass, which depends on the conversion mechanism. 
Let us now look at the other observables. 

Gravity waves. Relic gravitational waves will have a very blue spectrum since they are just 
sensitive to the ratio H/Mp during the contracting phase and will be therefore strongly suppressed 
on large scales. The situation is very close to ekpyrotic/cyclic models [23] and no detection is possible 
in the foreseeable future. 

Non-Gaussianities. There are two sources of non-Gaussianity one has to consider. First of 
all, the exponential potential is non-linear and therefore it sources interactions among different 5ip 
modes. This is similar to what happens in inflationary models with a second field: self-interactions 
of this isocurvature scalar induce a non-linear evolution, so that the statistics of the field becomes 
non-Gaussian. The level of non-Gaussianity, i.e. the correction to the linear theory, can be obtained 
comparing the mass term ^V" (ip)5ip 2 with the cubic term ^ V'" (^Si^ . This gives 

NG^^- (18) 
MM y ' 

As 5ip increases with time, most of the non-linearity will develop at the very end. The resulting 
non-Gaussianity is thus of order l/(|to|M), where to is the latest time where the solution described 
in the previous Section applies. 

The leading non-Gaussianity is cubic (quartic and higher terms are suppressed with respect to 
the cubic by further powers of 5ip/M) and it will show up as a 3-point function correlator. As 
non-linearities are dominated by the late time evolution, when gradient terms can be neglected, the 
momentum dependence of the 3-point function will be of the "local" form. The same happens in 
multi-field inflationary models where the out of the horizon evolution dominates the non-Gaussianity 

[22 ES]. 



3 In the text we consider the case in which the ratio H 2 /\V"(il>)\ is small and constant, which is analytically 
simpler. However it seems rather tuned to assume that the two quantities H 2 and V" vary by many orders 
of magnitude keeping their ratio constant. To have approximate scale invariance we just need this ratio to 
remain small, though it can vary significantly while modes of cosmological interest freeze. 
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The explicit calculation, neglecting 0(p) corrections, is quite simple. Following Maldacena [27] 
the 3-point correlator at the final time to is given at tree level b}0 

(6^0)6^0)6^0)) = -i [ ° (6^0)6^0)6^^0) H- mt (s))ds + c.c. (19) 
where the interaction Hamiltonian is the cubic 6ip self-interaction 



Using the normalized free field solution 



1 / ^ i 



<% = -^e~ lkt ( 1 - 7- ) (21) 



we get 



(%J(o)% J (to)%,(«o)) = (2^) 3 ,S(^fc,)=^ p (J- +ih) (i+ife) 



We are interested in the leading contribution for small to, which is given by 

We recognize in front the expected "local" momentum dependence. With the usual definition of 
the power spectrum of 6ip at time to- (6ip^ 6ip^) = (2ir) 3 6(ki + fo)^ 3 ■ A^, A^ = 1/(2^), we can 
rewrite the result as 

<%>)^>)% 3 (*°)> = ( 2 -) 3 ^(E^)^f A f • ' (24) 

We see explicitly that the non-linear corrections are of order l/(M|to|) as discussed above. 
The local form of non-Gaussianity is usually defined through the relation 

C(x) = Cg(x) - lfhT%(x) 2 ~ (Q) , (25) 

where £ is the observed perturbation and Q g is a Gaussian variable. Experimental limits are given 
in terms of the scalar variable f^i^- The variable Q will be related to Sip as £(af) = ±6ip(x)/M c , 
with M c a mass scale depending on the conversion mechanism. Notice that we have a sign ambigu- 
ity: depending on the mechanism, positive Sip will correspond to positive £ or viceversa. Possible 
additional non-linearities in this relationship will be discussed below. 



4 The tree level calculation can be interpreted as the classical non-linearity among the modes generated by 
the background. There will be quantum corrections to this calculation coming from loop diagrams and thus 
suppressed by higher powers of h [28]. In our case an interaction of the form V^ n 'Sip n /n\ will give a correction 
to the 3-point function of order 1/ (tM) 2 (E /M) n ~ A , where E is the typical energy of the process E ~ |*| — ^ 
The cubic term is the most important giving a correction ~ 1/ (|f |M), which is anyway small as the tree level 
calculation forces l/(|fo|M) <C 1. 
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The 3-point function of £ will thus be given by 
which implies 



f local _ 5^ 1/2 1 , . 

~ 24 ^ M\t \ • ( } 

The tightest experimental constraints on are presently coming from the analysis of WMAP 



3yr data [2U |2g] 
This constraint gives 



36 < /nl < 100 at 95% C.L. (28) 



1 24 i/o o 

— — r < 100 = Ay ~ 6 • 10~ 2 at 95% C.L. (29) 

M|t | 5^ C 



The second source of non-Gaussianity is the conversion between isocurvature perturbations 5ip 
and the adiabatic mode Q. As everything happens when gradients are irrelevant the relationship 
between the two variables will be local in real space, but in general not linear 

C(f) = 5^{x)/M c + f^\5^{x)/M c ) 2 + ... (30) 

The quadratic correction gives an additional contribution to the local non-Gaussianity calculated 
in eq. (|27p . What is experimentally constrained is the sum of the two parameters f^i^- The 
contribution from the conversion mechanism cannot be further studied without specifying the model; 
only the local shape dependence of the 3-point function is model independent. 

In this Section we have derived all the predictions assuming that the energy density of the field ip 
is so small that its mixing with gravity can be neglected, at least when modes relevant for cosmology 
freeze. In Appendix Owe study the full action describing both ip and the field cj>, which dominates 
the stress energy tensor. In this way we can check under which conditions the mixing of tjj with 
gravity can be neglected. If the mixing cannot be disregarded, predictions will change and become 
more model dependent, as they will depend on the Lagrangian of <f). In this case our sharp prediction 
of a blue tilt might also change. We leave these scenarios for future study. 



3 The whole story 

Now that we have discussed the dynamics of perturbations sourced by the ip field, let us concentrate 
on the unperturbed history of our bouncing cosmology. We will assume that the entire history, until 
the beginning of the standard hot FRW era, is dominated by the scalar cj>, with the Lagrangian of 
a ghost condensate theory with the addition of a suitable potential. As it should be clear from the 
discussion above, this model should be regarded as an explicit example of a smooth evolution from 
contraction to the bounce and then to the expanding phase, which is not affected by fast instabilities 
or other pathologies. As predictions do not depend on the explicit realization of the unperturbed 
history, alternative models may be found, either at the level of field theory or in which stringy effects 
are important. 
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3.1 Contracting phase 

Although not required, we saw in the previous Section that it is simpler to study an evolution with 
a ~ \t\ p with p constant. Furthermore we demand p <C 1, so that the spectrum is close to scale 
invariance. We want to find out which form of the potential V(4>) must be added to the ghost 
condensate Lagrangian to obtain this evolution of the scale factor. Neglecting terms with more than 
one derivative acting on (ft we have 

C = ^~g [P(X, (ft) - V(</>)] X = -g^d^dA ■ (31) 

Notice that we are discussing a slightly more general situation with respect to the Introduction, in 
which the shift symmetry of (ft is not only broken by the potential, but also in the generalized kinetic 
term P. The reason for this will become clear later on. For concretness we can take a functional 
form P(X, (ft) = (ft~ 2 P(X). In the absence of a potential, it is easy to check that, if P has a minimum 
where it vanishes (P'((fto) = P(4>q) = 0), we have the solution (ft(t) = (ftot with a Minkowski metric 
(or de Sitter in the presence of a positive cosmological constant). In this case the system behaves 
as the original ghost condensate, with the only difference that, as P t xx(<fto, <fto(t)) = (<ftot)~ 2 P" ((fto) 
depends on time, the dynamics of perturbations progressively changes. 
The stress energy tensor is given by 

Pjiv = [P(X, (ft) - V{(ft)]g^ + 2P(X, 0), x d^d u (ft . (32) 

For sufficiently small departure from P' = 0, the part of the stress energy tensor that depends on 
the generalized kinetic term is linear in 7r and it contributes only to the energy density, without 
pressure. In the Friedmann equation for H the contribution of the potential cancels: 

1 11 M 2 - . 

M 2 H = --(p + p) = -- • 4P XX (ft 4 ir = "2-^ , K = 4 ^>o • (33) 

In the last equality we have assumed that the deviation from the solution in the absence of potential, 
(ft{t) = (ftot, is small. Assuming this, we can solve for tt as a function of t: 

- Ml^r = f^TT 7T = tL p . (34) 

F t 2 2 t 2 M 2 

The velocity of (ft is shifted by a constant amount. The deviation from the solution in the absence 
of a potential is therefore small if -k = 2p ~ 2 P < 1, as we are going to assume in the following. We 

gc 

can now solve the first Friedmann equation to find V((ft) 

V(t) = 3M 2 H 2 + 2M 2 H = M p( 3 ^ ~ 2 P) . (35) 
In the approximation tt <C 1 we have <p(t) = (ftot so that the potential turns out to be 

y __(»-V)^8 +ow . (36) 

(ft z 

We are interested in the regime p <C 1 so that we obtain an inverse quadratic potential which is 
negative and going to — oo for t — ► 0. This solution makes sense. In the absence of a potential, 
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starting from tt > 0, w6 would have a matter-domination like solution, as the contribution to 
proportional to tt does not have pressure. This would correspond to p = | and in fact V = in the 
equation above for p = 2/3. As t~ 2/ fr oc H oc a _2 / p , for smaller p we need t _2 7r to increase faster as 
a function of the scale factor and this corresponds to a negative tilted potential which induces an 
increase in the 4> velocity. Viceversa for p > | . 

It is important to stress that the solution above is a dynamical attractor in the same way as in 
the ekpyrotic/cyclic models. One can check that if one starts with a 7T slightly different from the one 
of eq. (|34p , this perturbation to the velocity dies as t and the solution goes back to the unperturbed 
one. Also curvature and anisotropies become irrelevant, as the energy density increases with the 
scale factor faster than a -6 [29]. For a local observer the evolution converges to a single unperturbed 
history, which turns out to be a crucial simplification to allow following perturbations across the 
bounce [TT]. This last point is not so relevant in our case as we have an explicit bounce and we can 
follow the perturbations throughout. 

To conclude this Section we are now going to study perturbations of <p\ the purpose is just to 
check that it is consistent to neglect them and to have only the i/j sector as source of the observed 
perturbations. 

Following the general results of [16J in Appendix [A] we derive the Lagrangian for the usual 
variable ( which describes the scalar perturbations of the system^]. We simply get 

S- J 2"^ C "if U C J 



(37) 



One can show by symmetry arguments pjj] that the coefficient of the spatial kinetic term (<9£/a) 2 is 
always given by MpH /H 2 , independently of the <fi Lagrangian. The time kinetic term is not fixed 
by symmetries and in our case it is time- independent; therefore £ modes have a constant speed of 
sound c s = v^Mp/Mgc <C 1. The amplitude of a £ mode at freezing will be 

C ^'S^^' (38> 

where t and H must be evaluated at horizon crossing. The amplitude grows as H so that the 
spectrum is strongly blue. Neglecting ^-corrections -freezing (&) oc k, so that the final spectrum goes 
as l/k instead of the usual k~ 3 for the scale invariant case. This also tells us that the C perturbations 
are completely negligible on scales of cosmological interest. 

As we discussed, in the limit H — > the spatial kinetic term in (|37p vanishes. For this reason it 
is important to consider higher derivative terms which in the ir Lagrangian give \17\ [T6] 

1 

~ 2 J 

This term is the leading spatial kinetic term in the original ghost condensate model for which H = 
and it is very important to control the stability of the system as we will discuss in the next Section. 

However, as a consequence of the mixing with gravity, the higher derivative term proportional 
to M 2 C induces an additional 2-derivative term in eq. (|37|) of the form 

M3 c M gc V 2 (0 y (40) 



-M 2 c (d\/a 2 ) 2 . (39) 



5 What is relevant for observations is the constant mode of ( [TT]. As we are interested in studying the 
behaviour of <f> perturbations, we simply disregard the tp field in this Section. 
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which has the wrong, unstable sign and describes a Jeans-like instability [T7]. The system is stable 
if the term (<9£) 2 in eq. (|37p . which has the healthy sign, dominates. This happens for 

M gc <C c s M P . (41) 

Restricting to this regime, one can also see that around the time when a mode freezes the term 
(d 2 X I a 2 ) 2 is completely negligible, and this justifies our use of the simplified Lagrangian ([37j) for the 
study of the fluctuations generated during the contracting phase. 



3.2 The bounce 

The very same procedure can be used to get the potential for the bouncing phase. For simplicity 
we come back to the usual ghost condensate with a non standard kinetic term of the form P{X), 
with the shift symmetry broken only by a potential term. In this case we can assume without loss 
of generality that AP"^ = <p 2 = M 4 C . 

As an example we take H to be a linear function of time across the bounce 

H(t) = A = T2- ( 42 ) 

As H is constant eq. ([33]) implies a constant negative tt so that 

2Ml . . 

v = r-^r* > (43) 

where we choose the integration constant such that tt vanishes when the Universe reverses from 
contraction to expansion. In the limit \tt\ <C 1, we get the parabolic potential 

2MJ 3Mj^ 
W) T 2 + (M gc T) 4 1 ' 

We start in contraction with a positive potential energy and a negative contribution from the tt. 
The potential decreases in time, while tt stays constant until the total energy density goes to zero 
and the Universe bounces to an expanding phase. 

The dispersion relation for the fluctuations of the ghost condensate as calculated in |16j for a 
generic FRW background is 

MlM^ + AM^H M 2 C 

As discussed in the last Section, there is a k 4 term proportional to M 2 C , and two terms going as 
k 2 . The difference with the last Section is that now both the k 2 terms have the unstable sign. The 
one proportional to M 2 C describes the Jeans instability which is a consequence of the mixing with 
gravity. The second term is proportional to H, as we remarked in the last Section, so that it has 
now the unstable sign. 

It is straightforward to check [16] that one cannot make both these sources of instability ar- 
bitrarily small: the best compromise, i.e. the choice which minimizes the instability rate of the 
system, is to make the two instabilities comparable and this happens for M gc M 2 c /Mp ~ H 1 / 2 . In 
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this case it is easy to check that the most unstable mode has a rate of instability ui max — 

As the duration of the bounce is of order T, this mode grows by order one and this is clearly not 

problematic. 

One may worry about the relevance of this instability for predictions, as modes of cosmological 
interest evolve during the bouncing phase. Intuitively however one expects that if the bounce occurs 
sufficiently fast, the time during which the very long wavelength modes evolve is much shorter than 
their typical frequency, making the effect of this evolution irrelevant for cosmological observations. 
Let us check that this is indeed the case and the effect is completely negligible. As one can easily 
verify, the most unstable mode has a wavelength much smaller than i/ _1 (T), so that for modes 
of cosmological interest, which are obviously with wavelengths much larger than i7 _1 (T), one can 
neglect the A; 4 term and take the dispersion relation to be 

Each mode will grow by a factor uj(k) ■ T, where we have taken the duration of the bounce to be of 
order T . Assuming that just after the end of the bounce we begin a hot FRW cosmology, a given 
mode can be characterized as in inflationary cosmology with its number of e-folds. In this context 
this just tells us how much its wavelength is larger than the horizon at the beginning of the FRW 
era. We obtain 

Mf c T M P v ; 

where in the second passage we have taken the relation T -1 = H 1 / 2 = M gc M 2 c /Mp, the choice which 
minimizes the instability rate. The effect is in any case exponentially small and thus completely 
irrelevant for observations. 

The two phases described above, the contraction with constant p and the bounce, must be 
somehow glued together. One can imagine for concreteness a sudden transition, when the potential 
for (p jumps from the negative values required for the contracting phase to the positive ones required 
for the bounce. In the transition, energy is conserved so that the raise in the potential must be 
compensated by a negative contribution: tt becomes negative, which is indeed necessary for the 
bounce to happen. A schematic representation of the potential is shown in Fig. [TJ At represents the 
interval of time during which the transition from the contracting to the bouncing solution occurs. 
Notice that one can make At much shorter than H , so that it can be treated as instantaneous for 
the cosmological history, and at the same time much longer than M _1 , so that the transition can 
be described within the regime of the effective theory. 

The sharp jump in the potential is rather unpleasant. However it is important to stress that 
something so abrupt is unavoidable in a bouncing cosmology where H must evolve in a short time (of 
order H -1 ) from a large negative value before the bounce to a large positive value in the expanding 
phase. 

Finally, as shown in Fig. [TJ we imagine that at the end of the bounce the energy associated with 
the ghost condensate is converted to radiation, and the standard epoch of cosmology begins. 

We have not yet discussed the conversion of 5ip perturbations to adiabatic ones. One can envisage 
many possibilities. We assumed that when relevant modes freeze, the field ip gives a sufficiently small 
contribution to the energy density to neglect its mixing with gravity (see Appendix |Aj). If this ceases 
to be true later on, a mixing between 5ip and £ will occur. Alternatively a mixing between (j) and ip 
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Figure 1: Schematic representation of the potential of 4> during the contracting and the expanding 



in the potential can source the conversion. In both cases an approximately scale invariant spectrum 
for (" is induced in a way similar to what shown in Appendix [A] and this will be preserved until 
the expanding phase: the general approach advocated in [11] can in fact be applied and explicitly 
checked in our smooth bouncing solution. Another possibility is that the bouncing phase itself 
depends on the isocurvature scalar tp, a sort of "variable bounce" in analogy with the variable decay 
scenario in inflation [31]. In this particular setting, one could imagine a connection between the 
sudden transition in the potential and some coupling between 4> and tjj. A further possibility is 
that 5ifj fluctuations are converted only in the expanding phase, as in the curvaton [30] and variable 
decay models [31]. Whatever is the mechanism, to find the final value of £ one has just to solve 
the unperturbed history for different values of ifi in the commonly used "parallel Universes" or 5N 
approach. This simplication follows from the fact that 5ip is relevant only when relevant modes are 
out of the horizon. 



In inflation the origin of an approximately scale invariant spectrum of perturbations can be traced 
back to the fact that the inflaton potential is very flat, which implies that the Hubble parameter 
is almost constant. In our scenario the origin of scale invariance is completely different, as it 
comes from a precise functional form of the potential V{ip), i.e. a negative exponential. It is worth 
emphasizing that while relevant perturbations are produced, the potential of ip changes by many 
orders of magnitude; an exponential is not a useful approximation to simplify the algebra but a 
physical request: we need a negative exponential potential over a large region. 

One could object that this is not very satisfactory: in inflation we just need a sufficiently shallow 
potential, while here we have to require a precise functional form, which includes an infinite series 



phase. 



4 Naturalness of the exponential 
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of non-renormalizable operators! We are obviously interested in the perturbative region of the 
exponential, where the potential flattens out and the theory becomes weakly coupled. The flatness of 
the potential in this region is technically natural, in the sense that loop corrections will be suppressed 
by derivatives of the potential, which are small. In other words ip can be thought as an approximate 
Goldstone boson in this region. However this does not help in selecting a specific function, among 
the infinitely many which flatten out at infinity. 

The point is that an exponential potential is not a generic infinite series of non-renormalizable 
operators, but a particular one which emerges in many examples. There are at least two possible 
origins of exponential potentials. The first one is just a consequence of dimensional reduction as it 
can be seen in this very simple example. If we take a 5-d theory compactified on S 1 , the dimensional 
reduction of the Einstein-Hilbert action gives at the level of zero modes 

S = \ f afx^MlTR, (48) 



2 _ 

where g is the 4d metric and R its Ricci scalar. The radion field T multiplies the 5d Planck scale 
M5. If we do a conformal rescaling to go to Einstein frame the radion gets a kinetic term of the 
form 

d 4 xV=gMl^ (dlogT) 2 , (49) 

where M4 is the 4d Planck mass. This implies that if the energy of the system depends polynomially 
on the radius T, for example including a cosmological constant or just at loop level because of 
the Casimir effect, the canonically normalized radion will have an exponential potential, with a 
positive or negative overall sign depending on the model. For instance a negative cosmological 
constant induce a negative exponential potential. The generalization of this trivial example leads to 
the appearence of exponential potentials in supergravity compactifications and therefore in string 
theory. 

Another possible origin of exponential potentials in through non-perturbative effects. In su- 
persymmetric theories one can have directions in field space which are flat in perturbation theory. 
These can however be lifted by non-perturbative effects. For instance if the flat direction enters in 
the gauge coupling, instanton effects will give contribution going as powers of 



(50) 



We thus see another possible origin of an exponential potential. 

We stress that, although an exponential shape of the potential is rather ubiquitous in many 
scenarios, it is not clear whether the condition M <C Mp and the request of a negative potential can 
be fulfilled in string theory [6'2\ 155] or in some other UV complete theory. 



5 Conclusions 

Inflation is a tremendously compelling scenario. However to make progress it is crucial to have 
alternative models and to work out their predictions, leaving the final word to experiments. It 

6 The situation resembles the one of inflationary models with a variation of the inflaton larger than the 
Planck scale. Without additional symmetries, we need a functional fine tuning to keep the potential sufficiently 
flat, i.e. we have to control an infinite series of non-renormalizable operators. 
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might be dangerous to think that experiments are just measuring the parameters of the only theory 
we have. 

In this paper we described a bouncing cosmology, in which the bounce is smooth and without 
pathologies. This is indeed possible in a theory with a scalar with a modified kinetic term, like the 
ghost condensate model with the addition of a potential. 

What is new in this scenario with respect to previous bouncing models, like the pre Big-Bang 
scenario and the ekpyrotic/cyclic one, is that the physics of the bounce is explicit and under control. 
On the contrary, in the other models one must assume that the bounce occurs and its description 
must remain qualitative as it lies out of the regime of validity of effective field theory and of the 
present understanding of string theory. This is particularly important when one studies the evolution 
of perturbations across the bounce, which is crucial to assess the viability of a model. To this smooth 
model of a bounce we add the mechanism to produce density perturbations recently studied by 
Lehners, McFadden, Steinhardt and Turok [18] in the context of the ekpyrotic/cyclic model, that is 
an isocurvature scalar moving along a negative exponential potential. 

We therefore have an explicit and controllable model, in which predictions can be derived: no 
observable gravitational waves, a high (but model dependent) level of non-Gaussianity with a local 
shape and, most importantly, a slightly blue spectral index. This model will be ruled out if the 
present preliminary detection of a red tilt is confirmed by future data. 

The model is clearly not very compelling as the Lagrangian of the <j> field seems really ad hoc. 
This is the result of the fact that it is difficult to construct a non-pathological system which induces 
a bounce. The stress energy tensor must violate the null energy condition and this usually leads to 
catasthrophic instabilities [16]. Maybe there are much simpler systems leading to a bounce, either 
at the level of effective field theory or in which quantum gravity is relevant. Notice however that 
predictions depend only on the isocurvature scalar tp and therefore are not sensitive to the explicit 
realization of the contracting phase and of the bounce. 

Alternatively, it might be that the complicated Lagrangian for (j) signals the fact that theories 
which violate the null energy condition cannot be realized, in the sense that our effective field theory 
cannot be UV completed. A deeper understanding of the implications of the null energy condition 
is clearly needed. 
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Appendix 



A Gravitational mixing between </> and ip 

In this appendix we provide the explicit Lagrangian for the system of the Ghost Condensate (j) plus 
the scalar field ip during the contracting phase. With the resulting Lagrangian, we shall be able to 
study the conditions under which we can neglect the gravitational fluctuations in the 5ip equation 
of motion, justifying more regourosly what assumed in sec. [2j 

In order to construct the action, we follow closely [27]. We use ADM parametrization of the 
metric: 

ds 2 = -N 2 dt 2 + g ij (dx i + N i dt)(dx j + N j dt) , (51) 
where our background solution is of the form: 

ds 2 = -dt 2 + a 2 (t)dx 2 , (52) 

where a(t) oc \t\ p , and we consider general metric fluctuations 

N = 1 + 5N , Nj = 5Nj , gij = a 2 (t)5ij + dcjij . (53) 

In this language, the Einstein-Hilbert action takes the form: 

Seu = \m 2 p jd A x^g-R= hi 2 J d 3 x dt sfg [NR® + -^(E* ' E {j - E\ 2 )} , (54) 

where is the Ricci scalar of the induced 3D metric. E^ is related to the extrinsic curvature Kij 
of hypersurfaces of constant t, 

E id = NKij = ^[d tgij - ViNj - V^Ni] , (55) 

where V is the covariant derivative associated to the induced 3D metric gtj. 

The ghost condensate <ft and the scalar field ip has a background solution of the form: 

4>{x, t) = 4> Q ■ (t + 7T (t) + Sn(x, t)) , (56) 

rP(x,t) = ijj (t) + Sip(x,t) . (57) 
The action for the Ghost Condensate is given b}Q: 

Sec = jd^xV^g P(X, <f>) = J d 3 x dt y^N P(X, <P) , (58) 



7 In this Appendix we neglect effects due to higher derivative terms that in the main text are proportional 
to Af g 2 c . Following [Tg], one can show that in order for a comoving mode to be stable before freezing 
(oj(k) ~ H), M gc must be so small that the k term becomes irrelevant well before this time. Furthermore, 
as discussed in the main text, the Jeans instability that M gc would induce is subdominant with respect to 
the k 2 term proportional to H . Here we are interested in studying the gravitational mixing between <f> and 
t/j, and this occurs only around or after the freezing time, so that we can safely neglect the higher derivative 
terms in our discussion. 
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where X = —g^ u d^pd v p. The action for the field ip is given by: 
1 f rr r „_i / : „ A 2 



S f = 2j d * X dt " iVia ^) " Ngf'diipdjij) - 2NVtyj\ , (59) 

where V(?/>) = Voe^/ M . We see that ?/> and (f> mix only through gravity. As we said in the main part 
of the paper, at a certain point during the contracting phase, or after the bounce, there must be some 
direct interaction between the <p and ip to allow for the conversion of the isocurvature fluctuations 
into adiabatic one. This however must occur at least after all the relevant cosmological modes have 
exit the horizon. This is a much later time with respect to the one we are considering here. For this 
reason, here we can neglect all direct interactions between ip and p. 

The ADM formalism is designed so that one can think of gij, p and ip as dynamical variables, 
and ./V and N l as Lagrange multipliers. We will choose a gauge for gij,4> and ip that will fix time 
and spatial reparametrizations. We find convenient to define the gauge: 

5vr = 0, gij = a 2 {t)(l + 2Q5 ij . (60) 

This gauge fixes completely time and space diffeomorphisms at non zero momentum. It represent a 
gauge where the ghost condensate is uniform, and taken as the time variable. The dynamical degrees 
of freedom are £ and 5ip. In the limit where the gravitational mixing of ip with <p is negligible, we 
can think of Q as the scalar degree of freedom associated with the ghost condensate. In this gauge, 
the action for <p takes the form: 

2 



S GC = I d 3 x dt 



-2M 2 + M 2 pP ) 1 t (l-3p) (-2M 2 +M 2 p) at | 1 M gc 



+ - -^-N + Ti-g-VN)' (61) 



t 2 N t 2 2 t 2 

The tadpole terms in (|6ip are chosen in order to ensure that the background solution has exactly 
the form a(t) cx \t\ p , and ipo(t) = — V(ipo(t)) = — 2M ^~ 3p ^ as the one we used in the main part 
of the paper. This means that the form of the function P(X, <p) used here and the corresponding 
4>o(t) solution will be slightly different (by order p corrections) from the one used in the main part. 
Thought the physical implications do not change, considering an exact solution simplifies relevantly 
the study of the action of the system. 

To find the action for £ and 5ip, we solve for iV and N l through their equations of motion and plug 
the result back in the action. This procedure is allowed because iV and N l are Lagrange multipliers. 
The equations of motion for iV and N l read at first order: 

\2M „ „ .i 

di [—^ + 2 M 2 P H{t) 5N-2 Mpq = , (62) 

2 Mpp(l - 3p) + Mj m 2M(1 - 3„) 2M 2M> _ _ I? = 

t 2 t 2 t t t a 2 

where for convenience we have defined: 

N* = diip , d 2 ip = x- (63) 
We can solve these equations to first order to obtain: 

Mpp p 



MM 2 M ■ 1 M 2 \ . t , o 



2M p p 2 t Y M 2 p Y \p 2M 2 p 2 j" pa 
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In order to find the quadratic action for ( and Sip, we do not need the second order solutions for 
N and N l . The reason for this is that the second order term in N will multiply the constraint 
^ evaluated at zeroth order, which vanishes since the zeroth order solution obeys the equations of 
motion. A similar argument holds for N % . 

Substituting in the action S = Seh + Sec + Ss-ip, after performing some integration by parts, we 
find an action of the form: 

S = Jd 3 x dt a 3 [C c + C^ + £ M ixing] , (65) 



where: 



£{= (2M|p + M gc ). 2 _M|M c| m 



2p 2 ' p \ a 

i ,-,2 i /a x 2 




MM„ C 2 ■ M ■ ■ M P fd 2 \ 

We see that, neglecting the mixing, the speed of sound of ( is given by: 

o 2pMl 2pM% 

c 2 = J- p , ~ JL-^L < l (67) 

{2Mlp + MD M 2 C 

where in the second and third passage we have used the fact that for the validity of the ghost 
condensate effective theory we need to have: 



pM 



M 2 



< 1 . (68) 



The speed of sound of C must be very small. It is useful to write the former Lagrangians in terms 
of c s . Keeping only leading order terms in c s , we obtain: 

^=H(c 2 -<!<) 2 ), (69) 



1-2 1M V / 6M 2 M 2 \ 1 2 

^ = 2^ "2U V V 1_3p "W + ' 

^Mixing = -2— 2-C<ty + 2 — + 2-f -^C W ■ 

If we look at the Lagrangian for #0, we see that the inclusion of gravitational perturbations has 
produced two effects: the appearance of a direct coupling between £ and tp, and a correction to the 
mass term which is generated when we plug back into the solutions for N and N l . Obviously, 
the form of the mixing terms does depend on the nature of the field which leads the contraction. 
However, there is a regime where ip is so subdominant that these mixing effects are negligible. In 
this case, the nature of the field which drives the contraction is irrelevant for most of the results. 
This is the case in which we concentrated in the paper. Now we are able to explicitly check what is 
the condition under which this is verified. The condition will turn out to be not exactly p^/ptot *C 1, 
because we shall have to impose the mixing to be so small not to alter the calculation we did for 
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the tilt, and also because the unmixed £ has a small speed of sound. It is rather straightforward, 
comparing for example the terms of mixing with the diagonal terms for 5ip, that the effects of 
gravitational mixing are negligible for the 5vp equation of motion if 

— = nSr « Mm{p 2 c;\pc 2 s } (70) 

It is interesting to note that even though this condition is satisfied, which means that £ has no effect 
on Sip, still at late times, t 3> 1/fc, the effect of Sip on Q is not negligible. This is just the effect that 
on large scales, after the particular mode has freezed, different values of Sip induce small differences 
in the expansion of the separate regions. The induced £ is scale invariant and of size: 

M l M 2 p 2 c 2 _ 5 

Cinduced~^< M ^« log < 10 , (71) 



where the first inequality comes from the limit on non-gaussinities, the second from (|70p . and the 
third from the fact that, because of the constraint on the tilt of the two point function, p < fewx 1CP 2 , 
and c s Cl for the validity of the effective field theory. The induced £ is therefore too small, justifying 
our request for a separate mechanism for conversion of isocurvature perturbations into adiabatic ones. 

An alternative interesting regime, different from the one on which we concentrated in this paper, 
is the one where the inequality (|70p begins to be violated. In this case, two important things might 
occur. On the one hand, it is clear from the mass term in that the factors of p that determine 
the tilt might be overcome, inducing possibly a red tilt. On the other hand, the induced £ might 
become large enough, without loosing its scale invariance, so that there might be no need for an 
explicit conversion mechanism between isocurvature and adiabatic perturbations. Though these 
results would be appealing, the predictions would depend also on the details of the field that drives 
the contracting phase, and therefore they would be much more model dependent. For this reason, 
in the main part of the paper we concentrated in the region where the condition (|70p is satisfied and 
all the effects coming from mixing with gravity can be safely neglected. 
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